A Bayesian approach to multiscale inverse problems using the sequential Monte Carlo method
نویسندگان
چکیده
A new Bayesian computational approach is developed to estimate spatially varying parameters. The sparse grid collocation method is adopted to parameterize the spatial field. Based on a hierarchically structured sparse grid, a multiscale representation of the spatial field is constructed. An adaptive refinement strategy is then used for computing the spatially varying parameter. A sequential Monte Carlo (SMC) sampler is used to explore the posterior distributions defined on multiple scales. The SMC sampling is directly parallelizable and is superior to conventional Markov chain Monte Carlo methods for multi-modal target distributions. The samples obtained at coarser levels of resolution are used to provide prior information for the estimation at finer levels. This Bayesian computational approach is rather general and applicable to various spatially varying parameter estimation problems. The method is demonstrated with the estimation of permeability in flows through porous media.
منابع مشابه
A Multiscale Strategy for Bayesian Inference Using Transport Maps | SIAM/ASA Journal on Uncertainty Quantification | Vol. 4, No. 1 | Society for Industrial and Applied Mathematics
In many inverse problems, model parameters cannot be precisely determined from observational data. Bayesian inference provides a mechanism for capturing the resulting parameter uncertainty, but typically at a high computational cost. This work introduces a multiscale decomposition that exploits conditional independence across scales, when present in certain classes of inverse problems, to decou...
متن کاملE-Bayesian Approach in A Shrinkage Estimation of Parameter of Inverse Rayleigh Distribution under General Entropy Loss Function
Whenever approximate and initial information about the unknown parameter of a distribution is available, the shrinkage estimation method can be used to estimate it. In this paper, first the $ E $-Bayesian estimation of the parameter of inverse Rayleigh distribution under the general entropy loss function is obtained. Then, the shrinkage estimate of the inverse Rayleigh distribution parameter i...
متن کاملBayesian analysis of the scatterometer wind retrieval inverse problems: some new approaches
The retrieval of wind vectors from satellite scatterometer observations is a non-linear inverse problem. A common approach to solving inverse problems is to adopt a Bayesian framework and to infer the posterior distribution of the parameters of interest given the observations by using a likelihood model relating the observations to the parameters, and a prior distribution over the parameters. W...
متن کاملSpatial count models on the number of unhealthy days in Tehran
Spatial count data is usually found in most sciences such as environmental science, meteorology, geology and medicine. Spatial generalized linear models based on poisson (poisson-lognormal spatial model) and binomial (binomial-logitnormal spatial model) distributions are often used to analyze discrete count data in which spatial correlation is observed. The likelihood function of these models i...
متن کاملBayesian analysis of the scatterometer wind retrieval inverse problem: some new approaches
The retrieval of wind vectors from satellite scatterometers is a non-linear inverse problem. A common approach to solving inverse problems is to adopt a Bayesian framework and infer the posterior distribution of the parameters of interest given the observations using a likelihood model relating the observations to the parameters, and a prior distribution over the parameters. In this paper we sh...
متن کامل